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Abstract. Measurement of Cosmic Microwave Background (CMB) anisotropies has been playing a 
lead role in precision cosmology by providing some of the tightest constrains on cosmological models 
and parameters. However, precision can only be meaningful when all major systematic effects are 
taken into account. Non-circular beams in CMB experiments can cause large systematic deviation in 
the angular power spectrum, not only by modifying the measurement at a given multipole, but also 
introducing coupling between different multipoles through a deterministic bias matrix. Here we add 
a mechanism for emulating the effect of a full bias matrix to the Planck likelihood code through 
the parameter estimation code SCoPE. We show that if the angular power spectrum was measured 
with a non-circular beam, the assumption of circular Gaussian beam or considering only the diagonal 
part of the bias matrix can lead to huge error in parameter estimation. We demonstrate that, at 
least for elliptical Gaussian beams, use of scalar beam window functions obtained via Monte Carlo 
simulations starting from a fiducial spectrum, as implemented in Planck analyses for example, leads 
to only few percent of sigma deviation of the best-fit parameters. However, we notice more significant 
differences in the posterior distributions for some of the parameters, which would in turn lead to 
incorrect errorbars. These differences can be reduced, so that the errorbars match within few percent, 
by adding an iterative reanalysis step, where the beam window function would be recomputed using 
the best-fit spectrum estimated in the first step. 
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1 Introduction 

Measurements of Cosmic Microwave Background (CMB) opened a whole new era in theoretical 
physics. Not only it made the standard cosmological model widely acceptable, but it offered pre¬ 
cise measurements of cosmological parameters. Several high precision ground based and space based 
experiments were carried out in the past decades to measure CMB anisotropies. Recently, WMAP 
[1-4] and Planck [5, 6] mapped the full CMB sky with few arcmin level resolutions. However, high 
precision measurements demand accurate accounting of the systematic errors. Systematic errors can 
arise at different stages in CMB analysis, such as, foreground cleaning [7, 8], instrument calibra¬ 
tion [9, 10], measurement of beam response [11, 12] and estimation of the angular power spectra [13] 
from the observed maps. In this paper we focus on the effect of noncircular beams [14-21], one of the 
most important and challenging sources of systematic error. 

Observed CMB sky is a convolution of the underlying true sky with the instrumental beam 
response function. Accounting for the beam is thus necessary for measuring the statistical properties 
of the true sky. This turns out to be trivial if the beam is symmetric about the pointing direction, an 
unbiased estimator of the true angular power spectrum can be readily obtained. However, an actual 
beam response function is never so symmetric. The intrinsic optics of the instruments, aberrations due 
to the placement of detectors away from the principle optic axis, non-uniform distribution of pointing 
directions in a pixel, finite sampling duration in scanning and many other effects can distort a beam 
to make it non-circular. The asymmetries become progressively important in different analyses as 
the experiments strive to extract almost all the information embedded in the anisotropies. Beam 
asymmetry may lead to spurious effects in the measured CMB data. It can bias the angular power 
spectrum in a non-trivial way by introducing coupling between different multipoles through a bias 
matrix [16]. It can also cause statistical isotropy violation in the CMB anisotropy maps [22-25]. 
Here we study how the distortions in the angular power spectra caused by non-circular beams affect 
cosmological parameter estimation. 

If the bias matrix that couples different multipoles was available, it appears obvious that one 
would use it to obtain an unbiased estimator of angular power spectra and perform parameter estima¬ 
tion starting from the unbiased estimator. However, this process would encounter two major hurdles. 
First, attempt to unbias the estimator will introduce covariances among different multipoles. Second, 
perhaps the more important one, it is highly non-trivial to obtain the bias matrix for an experiment 
with non-circular beam shapes and complex scanning strategy. In literature the near-diagonal com¬ 
ponents of the bias matrix has been computed to leading order for simple scanning strategies only for 
temperature anisotropies [16, 18], though efforts are on to extend the method to polarization also [26]. 
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Our present work thus, in fact, evaluates the need for investing (enormous) effort in getting the full 
bias matrix for experiments like Planck in order to perform unbiased analyses. 

Here we develop a mechanism to include the effect of a bias matrix in the parameter estimation 
code SCoPE [27] in conjunction with the Planck likelihood code [13, 28]. We use Planck likelihood 
code in order to get realistic representation of noise. However, instead of deconvolving the observed 
power spectrum to get an unbiased estimator, we apply “inverse distortion” to the theoretical spectra 
in the likelihood code to get the desired posterior distributions. This step alleviates the need for incor¬ 
porating a non-trivial covariance matrix. We use different scalar transfer functions to do parameter 
estimation from the same convolved spectra and study the differences in posterior distributions with 
the correct estimates. 

It is very common in CMB analysis to use an effective scalar window function, which inherently 
assumes the beam to be azimuthally symmetric about the pointing direction. However, if the effective 
window function was derived using Monte Carlo simulations [19], including the full details of beams 
and scanning strategy, as was done for Planck analyses, one may be able to deceive an analysis by 
closely emulating the effect of a non-circular beam, as long as the fiducial power spectrum used for 
the simulations is close to the true one. Here we use WMAP and Planck best-fit power spectra to get 
two different scalar transfer functions and use them for estimating parameters from a test observed 
power spectrum obtained by convolving the Planck best fit spectrum with a bias matrix. However, 
since a full bias matrix for a mission like Planck does not exist, for numerical computation we must 
limit ourselves to a case where it is available. We use an elliptical Gaussian beam of similar size and 
ellipticity as one of the Planck high frequency detectors and non-rotating scan pattern. 

This paper is organised as follows. We present a brief primer on the connection between non¬ 
circular beams and window function in section 2. The general strategy for studying the effect of 
non-circular beam on cosmological parameters and the results of our study are included in section 3. 
Section 4 presents the conclusion and discussions. 

2 Effect of instrumental beams on CMB angular power spectrum 

Different experiments observe the CMB anisotropy field by scanning the sky through an instrumental 
beam of finite resolution. The observed time ordered data (TOD) is passed through a refined pixel 
binning procedure which is generally kept unaware of the beam shapes. This is because deconvolution 
of all the TOD samples (few trillions for Planck detectors) is computationally prohibitive. The 
resultant map can then be expressed as a convolution of the true sky with an effective beam function. 
If the effective beam for every direction was the same and symmetric about each direction, one could 
show that the angular power spectrum of the measured map Ci is trivially biased, Ci = Bf Cj, where 
Ci is the angular power spectrum of the true sky and Bi is the Legendre transform of the azimuthally 
symmetric beam function. It is indeed very common in CMB analysis to use an effective scalar window 
function Bf. However, perhaps for every CMB experiments, beams are non-circular. Non-circular 
beams make the effective window function tensorial, Ci = Aw Cy , as we briefly review below. 

2.1 Observed angular power spectrum 

The observed CMB anisotropy AT(q) in a given direction q can be expressed as the convolution of 
the true sky AT(q) with the instrumental beam function H(q,q') plus additive noise n(q), 

AT(q) = [ dOq' .B(q, q') AT(q') + n(q) , (2.1) 

J s 2 

where clflq is an infinitesimal solid angle around the direction unit vector q. Using the Spherical 
Harmonic transform of the observed map 

aim ■= I dOq Yj m (q) AT(q) = [ dflq [ dUq/Y) m (q)[H(q, q') AT(q') + n(q)] , (2.2) 

J s 2 J s 2 J s 2 
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and assuming Statistical Isotropy (SI) of the true sky, it can be shown [16] that the expected power 
spectrum of the observed map can be expressed as 

(Q) = 2ZTT ^ /, f (^(qi) AT(4a)) fl(qi ■ fc) (2.3) 

m=—l S 2 

= Y^AtvCv + C ; N . (2.4) 

l '=o 

Here, l max is the maximum harmonic multipole, C ; N is the angular power spectrum of noise and the 
bias matrix, 

A W ■= j d ^ 4l J dHq 2 P/(qi • q 2 )^'(qi,q 2 ) , (2.5) 

where, 

Wi(qi,q 2 ) := [ clHq/ [ dH^H(qi,qi)H(q 2 ,q2)P;(qi,q2). (2.6) 

J s 2 J s 2 

It is interesting to note that the noise-free two-point correlation function of the observed CMB 
anisotropy sky can be expressed as 

<AT(q)AT(q')> = £ (2 ^Q IH ; (q, q'). (2.7) 

;=o 


2.2 Circular Beams 

It is fairly straightforward to show that if the beam is azimuthally symmetric about the pointing 
direction, that is, i?(q, q') = B( q ■ q'), so that the beam can be expanded in terms of Legendre 
polynomials, 

^max 

B{ q,q') = B( q • q') = — ^(2Z + 1) B, fl(q • q') , (2.8) 

n 1=0 

the observed angular power spectrum is trivially biased, Aw — SwBf , 

Ci = BfCi + Q N . (2.9) 

Thus in this case it is easy to get an unbiased estimator using the scalar window function Bf , 

Cfcalar = B~ 2 [Ci - Cf] , (2.10) 

assuming that the noise power spectrum, C) N , can be precisely estimated independently from instru¬ 
ment noise characteristics. The co-variance of the unbiased estimators is given by 

Cov(Cf calar , ( 7 f/ calar ) = + B~ 2 Cf) 2 . (2.11) 

The above equations imply that in case of circular beams, there is no coupling between power spectrum 
at different multipoles. 

Note that, unless otherwise mentioned, we will use a twiddle (~) to denote an observed quantity, 
a bar ( “ ) on a quantity to denote an estimator and a superscript “scalar” to denote that the estimator 
is derived using a scalar transfer function. 

2.3 Non-circular Beams 

In general, however, B(q,q') can not be expanded in terms of Legendre polynomials, they can be 
expanded in terms of Spherical Harmonics Yj m (q') for every pointing direction q. Making use of 
Wigner rotation matrices, the spherical harmonic transforms of the beams computed at a specific 
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location could be transformed to other directions [14]. The final expression for the observed power 
spectrum can then be expressed in terms of the spherical harmonic transforms of the beam (6; m ) 
pointing at the 2 -axis of spherical polar coordinates. Numerical computation then shows that, for a 
trivial scanning strategy and elliptical Gaussian beams, to the leading order the bias matrix Aw has 
a large number of small off diagonal components, which can imply a significantly large difference in 
power spectrum if estimated assuming a circular beam [16]. Hence, from Eq. (2.4), one should define 
the true unbiased estimator as 

Cl = EAT' 1 ^' - c *]- (2- 12 ) 

l' 

Even though the approximate calculations provide a good estimate of the level of the effect, in 
order to fully account for the effect one needs the true bias matrix computation involving the actual 
non-circular beam shape, the exact scan pattern and beyond leading order computation. This is a 
challenging task and has not yet been accomplished in literature. Efforts are on to compute effective 
bi 2 for a complex experiment, incorporating the effect of scanning strategy, in turn enabling one to 
compute the true bias matrix to a reasonable accuracy using only leading order computation [29]. 

2.4 Scalar effective beam window functions 

Above we described that circular beam window functions can lead to significant systematic error in 
power spectrum estimation. But, realistic bias matrices for non-circular beams are not available. This 
situation might seem discouraging ! In practice though there is a middle-ground. One can compute a 
scalar beam window function through Monte Carlo simulations, incorporating beam asymmetry and 
scanning strategy. A large number of maps are simulated from a fiducial angular power spectrum, 
which are convolved with actual non-circular beams. The scalar window function can then be esti¬ 
mated by taking the average ratio of the convolved to unconvolved power spectra. However, a scalar 
window function corresponds to a circular beam (bi m oc <5 m0 ), as non-vanishing bi m lead to off-diagonal 
terms in the bias matrix. So this method essentially replaces the complex effective beams, which can 
also vary across the sky, by one effective circular beam. The goal of this paper is to study the accuracy 
and adequacy of this scalar beam window functions in cosmological parameter estimation. 

The effective beam window function is defined as the average of the ratio of the convolved (Ci) 
to unconvolved fiducial (C) fid ) simulated maps 

Wi := (a/Cf d ). (2.13) 

In this paper we simulate a convolved power spectrum by multiplying an unconvolved spectrum (Ci) 
with a given bias matrix 

Ci := Y^AwCv. (2.14) 

v =o 

This method alleviates the need for performing the Monte Carlo simulations for estimating the scalar 
transfer function, as neither Aw nor the fiducial spectrum Ci are random realizations. However, this 
is possible here because we are restricting ourselves to a case for which the bias matrix is available, 
which is a primary requirement for this study. 

Note that, in this paper we will be using two different unconvolved spectra, one is the “true” 
spectrum of the sky, which the analysis aims to recover, the other one is a “fiducial” spectrum (Cf d ) 
for estimating the scalar window function via Eq. (2.13). In both the cases convolution is done through 
Eq. (2.14). If one includes the the full bias matrix in the unbiased estimator, as in Eq. (2.12), the 
estimator is truly unbiased. This is because, combining with Eq. (2.4), one gets 

(Ci) = l(Cv) - C?] = Ci. (2.15) 

v 

However, this need not be the case if one uses a scalar transfer function Wi to define the unbiased 
estimator, as in Eq. (2.10), because 

<C? calar > = Wr'K Cl) - q N ] = W- 1 f I Aw Cl,. (2.16) 

J'=o 
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Inserting the convention for “true” and “fiducial” described above one gets 


(C7 Cal ar) = C\ 


fid 


Ei'=i 


max A 
0 


r true 
V 


Er=o All, cf, d 


(2.17) 


It is easy to see from the above equation that if Cf d = C) true , one would get ((7 ; scalar ) = C) true . However, 
this still does not ensure that the posterior distribution of the parameters can be correctly recovered 
by using the best-fit C) as Cf A . This is because parameter estimation codes compute likelihood for 
different sets of C) true , obtained by sampling the parameter space, to get the posterior distributions. 
Hence Cf d can match C) true at most for one realisation of the set of parameters. Moreover, since the 
aim of an experiment is to estimate the best-fit Ci one can not guess Cf rue a priori. Hence, there is no 
reason to assume that Cf d = C) true . Nevertheless, the difference between Cf d and C) true need not be 
very large either, as we do have approximate knowledge of C) true from previous experiments and only 
the “nearby” regions of the parameter space are sampled. In this work, as a realistic test case, we 
first take Cf d to be the WMAP best-fit spectrum and C) true to be the Planck best-fit spectrum and 
study whether this difference causes significant deviation in cosmological parameter estimation. We 
then study whether the true posterior distribution can be recovered (not only the best-fit parameter 
values) by taking Cf d and C ; true to be the same, the Planck best-fit power spectrum. 


3 Analysis and Results 

The broad approach we follow here is to generate non-circular beam convolved power spectrum, 
analyse it including the corresponding bias matrix in the parameter estimation codes and compare 
the results with the ones obtained using different scalar transfer function. 

3.1 Generation of Convolved Spectrum 

As mentioned above, we generate a convolved singular power spectrum Ci starting from a fiducial 
spectrum Cf d by multiplying the later with the bias matrix Aw [Eq. (2.14)] corresponding to the 
chosen non-circular beam. However, the bias matrix Aw has so far been computed to leading order 
for very specific cases of elliptical Gaussian beams. We use one such bias matrix for a beam of size and 
shape crudely similar to that of one of the Planck beams. Note that for this study it is reasonable 
to assume an approximate bias matrix to be the true bias matrix as the scalar window function is 
estimated through the same bias matrix. 

We choose the Full Width at Half Maximum (FWHM) of the beam to be 0.1° and eccentricity 
of e = 0.7. In literature often an elliptical Gaussian beam is characterized by ellipticity, the ratio of 
the semi-major (a) to the semi-minor (6) axis, while eccentricity is defined as e = -^/l — b 2 /a 2 . Hence 
the ellipticity of the beam used here is 1/v/l — e 2 = 1.4, which is on a slightly higher side compared 
to the Planck detectors, leading to a more conservative estimate which will become evident later. 

We chose a trivial scan pattern that keeps the alignment of the beam fixed with respect to the 
local meridian over the whole sky. We could not include a realistic scan strategy as the bias matrix for 
such a case has not been numerically computed in literature. In fact one of the main aim of this work 
is to study the need for investing enormous effort in computing the realistic bias matrix. However, this 
is not a big assumption either for this work. We are studying whether the scalar transfer functions 
can mimic a typical bias matrix in parameter estimation. Had we incorporated the full scan in the 
bias matrix, the matrix would be different, but we believe the broad characteristics would remain the 
same, as can be seen in the bias matrix plots for two different toy scan patterns presented in [16]. 
Hence this study should remain valid for realistic scan patterns. 

3.2 Choice of Scalar Transfer Functions 

We study four different variants of scalar transfer functions (IT)) for calculating the estimator of the 
power spectrum C) scalar ;= Ci/Wi as listed below: 
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1. Wi is taken as Bf, where Bi is the Legendre transform of a circular Gaussian beam, as defined 
in Eq. (2.8), whose radius is the geometric mean of the semi-major and semi-minor axis of the 
chosen elliptical Gaussian beam, preserving the total collection solid angle of the beam. 

2. Wi is taken as the diagonal components of the of the bias matrix (Aw), i.e. W) = An. 

3. Wi is computed from a “fiducial” Ci using Eq.(2.13). We use two different Cf d 

(a) Cf d is taken as the best fit WMAP-9 Ci, 

(b) Cf d is taken as the best fit Planck Ci. 

3.3 Modification to the Likelihood Analysis 

We compute the posterior distributions of cosmological parameters from the original Ci and the 
power spectrum estimators (C/) obtained by using different scalar and tensor transfer functions de¬ 
scribed above. Parameter estimation is done using the code SCoPE [27]. We add log-likelihoods 
from commander_v4.l_lm49.clik, lowlike_v222. clik and CAMspec_v6.2TN_2013_02_26.clik [28]. 
The effect of beam non-circularity is insignificantly small at the low multipoles. Since commander 
and lowlike only use C)’s from low multipoles, modifications to Ci with bias matrix or either of the 
above scalar transfer functions have negligible effect on these likelihoods. The only likelihood that 
gets affected in this process is that from CAMspec. 

CAMspec likelihood provides y 2 = (Ci — C z scalar ) T [C'] -1 (Cz — (7 ; scalar ) j where [C] is the covariance 
matrix and Ci is a theoretical spectrum for a given set of parameters. In most situations, to perform 
parameter estimation for a specific cosmological model, one runs MCMC chains over the parameter 
space and for each set of parameters C\ and y 2 are computed. This procedure does not question 
the validity of the estimator C) scalar and its covariance, which is in direct contrast to the situation 
considered in this work. Here we would like to modify C) scalar to capture the effect of using different 
scalar and tensor beam window functions. However, if we modify C'f calar by say W^ 1 A;z/Cf, calar , 
the covariance matrix also gets modified to [W i -1 A H /][C'][W z -1 Azz/] T . Therefore, for running MCMC 
with CAMSpec likelihood we need to change both C) scalar and the covariance matrix and feed those 
into the CAMSpec likelihood code. Rather, we use an alternate approach. For each MCMC realisation 
we multiply the Ci s obtained from CAMB[30] with Wi A^, , that is, replacing theoretical Ci by 
Wi ’Ey AJl'Ci' keeping C) scalar and [C] fixed. This step provides us the intended y 2 at a reduced 
complication. We then estimate the posterior distributions of cosmological parameters from the exact 
estimator Ci incorporating a bias matrix and the estimators C l i scalar for different approximations to 
the beam window described in Section 3.2. 

Note that, the above procedure is justified since the aim of this paper is not a reanalysis of 
Planck data, which would require a nearly exact treatment of beams and noise, the aim here is 
to verify if the scalar transfer functions are adequate for representing the distortions described by 
a bias matrix. So the above procedure in a way assumes that Planck measurements are correctly 
represented by the supplied beam transfer functions and noise covariance matrix and the default 
likelihood produces correct values. We distort the theoretical Ci by an inverse bias matrix and check 
if the effect can be compensated in the likelihood codes by multiplying it with a scalar transfer function. 
We use Planck likelihood code here only to get realistic noise characteristics and resolution. If one 
wanted to reanalyse Planck power spectra with bias matrices (assuming that they have somehow 
been made available for each channel), one would have to modify the observed spectra and the noise 
covariance matrices separately for each frequency. 

3.4 Results 

Wi \ ' ^—1 ^Planck 

We start by plotting the quantity Fi := — ' l ' - in figure (1), where C ; planck is the Planck 

best fit Ci and W) is estimated in three different (inexact) ways. This plot shows how much error 
is introduced in the power spectrum estimator for not considering the non-circularity of the beam 
properly. It can be seen that if we consider the beam to be circular Gaussian then the error in Cj 
at high multipoles is very large, more than ~ 3%. However, if we consider the diagonal components 
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of the Aw as the scalar transfer function then the error involved at high multipoles is lesser, ~ 1%. 
These levels of error in Ci can affect the cosmological parameters very significantly. On the other 
hand Wi computed through the forward method from the fiducial Ci, as in Eq.(2.17), performs much 
better. In the plot we choose the “fiducial” C; as C; WMAP . Here the error at high multipoles is less 
than 0.1%. Of course, if Wj was estimated using C' ; planck as the fiducial spectrum, T] would be 1 at 
every l. However, as mentioned before, C planck is used as the true spectrum here, which the analysis 
aims to recover, it can not be guessed a priori. We could also use some other set of parameters to 
generate the “true” spectrum. The Planck best fit parameters are chosen here as they promise to 
be the closest to reality in the history of CMB measurements. 



Figure 1. Comparison of different effective window functions. A convolved CMB power spectrum is 
generated by multiplying the Planck best-fit with a bias matrix for an elliptical Gaussian beam. It is then 
corrected by a scalar transfer function and its ratio to the original Ci is plotted. The scalar transfer function 
is obtained in different ways: by assuming a circular Gaussian beam (red), the diagonal of the bias matrix 
(blue) and using a forward approach with WMAP best-fit as the fiducial Ci (green). If we used Planck 
best-fit as the fiducial spectrum, the correction would be perfect and we would get Fi — 1. The good news is 
that even with WMAP best-fit fiducial Ci, Fi is close to unity. The other two scalar transfer functions lead to 
significant deviation of Fi from unity, implying highly incorrect estimation of the unbiased power spectrum. 

We have done parameter estimation for the standard six parameter LCDM cosmology, namely 
{fl c ,fl b ,h,T,n s ,A s } using SCoPE [27]. In figure (2) we show the likelihood contours (in blue) for 
the six parameters obtained from Planck best fit Ci, without introducing any instrumental effect. 
According to the scheme we followed to introduce the bias matrix, this case is equivalent of accounting 
for the full bias matrix in the analysis through the exact unbiased estimator [Eq. (2.12)]. Hence, this 
is the “correct” distribution we aim to recover. The contours corresponding to circular Gaussian beam 
approximation, the first C) scalar estimator discussed in Section 3.2, are overlaid (in red). It can be 
clearly seen that circular Gaussian approximation leads to very different posteriors, in some of the 
cases the contours get shifted by more than 2tr, implying that the mean values of the parameters 
can be highly incorrect as well. Using Wi = An also leads to similar, but little better, results. The 
likelihood contours are not shown here, the marginalised distributions are shown in figure 4 (in blue). 

Planck data analysis uses scalar transfer functions obtained via Monte Carlo simulations start¬ 
ing from a fiducial Ci as mentioned in Section 3.2. We first use the best fit WMAP Ci as the fiducial 
Ci and use the resulting scalar transfer function corrected power spectrum C) scalar for parameter esti¬ 
mation. We find that though the mean parameters are recovered to few percent of sigma accuracy, the 
distributions are significantly different, which would lead to wrong errorbars. Likelihood contours are 
not shown for this case, but marginalised distributions are plotted in figure 4. Next we show results 
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Figure 2. We plot the likelihood contours for the six cosmological parameters starting from the correct 
(unbiased using full bias matrix) pseudo-C; estimator Ci and a scalar transfer function corrected one (<5f calar ) 
obtained for circular Gaussian beam. It clearly shows that circular Gaussian approximation for a non-circular 
beam can lead to very large deviation in the posterior distribution. The correlation coefficients between pairs 
of parameters are shown in upper triangle of the plot. 


with Planck best-fit as the fiducial spectrum. In figure 3 the likelihood contours obtained from the 
full bias matrix corrected C) are plotted in blue and those from the estimator C) scalar are in red. The 
contours are close for most parameters, but not for all, the 1-D marginalised distributions provide a 
clearer picture. 

In figure 4 we overlay the 1-D marginalized probability distributions for the four different cases 
discussed in section 3, along with the distribution from the “correct” bias matrix corrected C; (green). 
It shows that if the WMAP best fit C) is used as the fiducial spectrum for deriving the transfer function 
(black) then also the average values of the distribution are almost the same, but the distributions are 
narrower, i.e. the standard deviations are smaller than that of the original distribution sometimes 
by as much as ~ 20%. This implies that the choice of wrong fiducial Ci for obtaining scalar transfer 
function may lead to wrong distribution of the parameter space, but not significantly changing the 
mean. When the fiducial spectrum matches the best-fit (yellow), the differences in distributions are 
reduced and the standard deviations differ by few percent, as in figure 3. 

Thus, we can conclude that the method adapted by the Planck team for obtaining the pa¬ 
rameters can provide reasonably accurate posteriors, at least if the beams are not too different from 
elliptical gaussian and the fiducial spectrum is close to the actual best-fit. However, the actual best-fit 
may not be known a priori. To address this issue, we suggest running the analysis first with a transfer 
function derived from a fiducial spectrum constructed from a reasonable initial guess of parameters, 
obtain the best-fit parameters, recompute the transfer function with the best-fit spectrum and finally 
get the posterior distribution and update the best-fit by using the new transfer function. 
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Figure 3. We plot the likelihood contours for the six cosmological parameters starting from the correct 
(bias matrix corrected) pseudo-C; estimator Ci (in blue) and a scalar transfer function corrected one (C'f calar ), 
where the transfer function was obtained using Planck best-fit as the fiducial Ci (in red). The correlation 
coefficients between pairs of parameters are shown in upper triangle of the plot. The plot shows that using a 
scalar transfer function derived via Monte Carlo simulations provides reasonably close best-fit values, which 
is clearer in the 1-D marginalized distributions. However, the posterior distributions do show deviations from 
the actual distributions, which become more significant if WMAP best-fit was used as the fiducial spectrum 
[see figure 4] instead of Planck , as the later may not be known a priori. These differences can be reduced 
by redoing the analysis using transfer functions derived from the best-fit spectrum obtained in the first step, 
the likelihood contours would then correspond to the red contours in this figure. 


4 Conclusions and Discussions 

We have estimated cosmological parameters from observed CMB temperature power spectrum by 
correctly accounting for non-circular beams with the help of a bias matrix and also by using different 
scalar beam window functions. We show that if the proper beam window function is not used, the 
power spectrum estimator gets biased and the posterior distribution of parameters becomes signifi¬ 
cantly different. Namely, assuming a circular Gaussian beam creates serious departure from the true 
posterior, while using a scalar transfer function estimated through Monte Carlo simulations involving 
non-circular beams does well in recovering the best fit values to few percent of sigma. However, 
the posterior distributions show significant differences, leading to incorrect errorbars, if the fiducial 
spectrum is very different from the actual best-fit. The posterior distributions can be brought closer 
to actual by introducing an iterative step, where the scalar transfer function is reestimated using the 
best-fit spectrum obtained in the first step. 
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Figure 4. One dimensional marginal probability distributions for all the four scalar transfer functions men¬ 
tioned in Section 3.2 are plotted along with the “correct/true” distribution that corresponds to the case where 
the bias matrix is accounted for. We show the average and the standard error below the plots. The deviation 
of the mean from the correct mean is shown in the bracket in terms of or, where a t is the standard deviation 
of the correct distribution for the parameter in context [that is, in the plot for n s , <tt is the standard deviation 
of the distribution for n s (green)]. Notice that both the yellow and the black marginalised distributions have 
means shifted only by few percent ctt from the correct analysis, illustrating that the scalar transfer functions 
derived through Monte Carlo simulations give accurate mean parameters, while the other approximations 
badly fail. However, the posterior distributions are different, the differences are more significant when the 
fiducial spectrum used for transfer function estimation is not close to the best fit (compare black and yellow 
with green). 


We have considered here an elliptical Gaussian beam and a non-rotating scanning strategy, as 
the bias matrix is available in only in such cases. If the bias matrix was available for a more realistic 
beam shape and scanning strategy, this work could be immediately repeated for that case. However, 
we believe that even in such complex situations, the bias matrix would be different but still would not 
show any dramatically different features. Since we have not used any special characteristics of the bias 
matrix, this work should still remain valid. Testing that hypothesis of course requires the whole bias 
matrix to be computed, though this work suggests that, at least for temperature and low-multipole 
polarization analysis, a full bias matrix may not be necessary for cosmological parameter estimation. 

Obtaining the bias matrix for polarization analysis is a even bigger challenge, though efforts are 
on [26] . We emphasise that the systematic effects of this kind become important when the errorbars 
are small even beyond the multipoles corresponding to the beam width. Which happens, for example, 
in CMB temperature anisotropy measurement with Planck. However, for polarization measurement 
the errorbars are still large, one may not need such fine corrections for current experiments. 

Here we would also like to mention that for our analysis we only change C ; TT , whereas a beam 
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function that affects Cf T should affect C EE . However the effect on Cj E is not significant in the 
present analysis, considering that the effects of the noncircular beam only affect the high multipoles 
and presently the polarization and the cross power spectrum measurement is limited low multipoles. 

It is well known that masking of foreground contaminated regions in the maps couples the lower 
multipoles. The full bias matrix for masking with circular beam is routinely used in CMB analysis [31]. 
Masking and non-circular beam effects are strong in different multipole ranges, still there can be small 
coupling between these two effects [18]. However, the effect of this on parameter estimation should be 
small and may not be significant at the level of precision attainable with current CMB experiments. 

The bottomline is that the need for incorporating full bias matrices in parameter estimation 
may not be crucial with current sensitivities of CMB measurements. Use of scalar transfer functions 
derived through Monte Carlo simulations, similar to what was done in Planck analysis, provides 
sufficient accuracy in best-fit parameter estimation and can also provide somewhat accurate posterior 
distributions if handled carefully. However, efforts should be invested in estimating the full bias 
matrices in these complex situations for putting more precise bounds on the cosmological parameters. 
Moreover, polarised bias matrices must be computed, to the leading order to start with, at least for 
performing a study like this one to make a statement about the adequacy of scalar transfer functions 
in high resolution polarisation analysis in Planck and post -Planck era. 
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